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Abstract —G.fast is an upcoming next generation DSL standard 
envisioned to use bandwidth up to 212 MHz. Far-end crosstalk 
(FEXT) at these frequencies greatly overcomes direct links. Its 
cancellation based on non-linear Tomlinson-Harashima Precod¬ 
ing (THP) proved to show significant advantage over standard 
linear precoding. This paper proposes a novel THP structure in 
which ordering of successive interference pre-cancellation can be 
optimized for downstream with non-cooperating receivers. The 
optimized scheme is compared to existing THP structure denoted 
as equal-rate THP which is widely adopted in wireless downlink. 
Structure and performance of both methods differ significantly 
favoring the proposed scheme. The ordering that maximizes 
the minimum rate (max-min fairness) for each tone of the 
discrete multi-tone modulation is the familiar V-BLAST ordering. 
However, V-BLAST does not lead to the global maximum when 
applied independently on each tone. The proposed novel Dynamic 
Ordering (DO) strategy takes into account asymmetric channel 
statistics to yield the highest minimum aggregated rate. 

I. Introduction 

Digital Subscriber Line (DSL) is dominant broadband ac¬ 
cess technology due to its ability to fulfill demands for 
reliable bigb-data-rate connectivity in a cost-effective way by 
exploiting tbe existing infrastructure of twisted-pair copper 
lines. Upcoming next generation DSL standard, G.fast H], 
proceeds in tbe trend of shortening copper lines (up to 250 m) 
between Central Office (CO) and Consumer Premised Equip¬ 
ment (CPE) aiming at fiber-like connection (up to 1 Gbps). 
Short lines enable the usage of wider bandwidth (initially up to 
106 MHz extended later to 212 MHz) than used by the forego¬ 
ing Very-high-bit-rate DSL (VDSL2) standard operating up to 
30 MHz. Cancellation of crosstalk between the lines by multi¬ 
user processing (denoted as vectoring or signal coordination) 
has a major impact on system performance. Ear-End-Crosstalk 
(FEXT) is the crosstalk affecting the other end of the line w.r.t. 
the transmitter as shown for Downstream (DS) in Fig.[T] FEXT 
is typically canceled by suitable transmitter precoding. If 
signal coordination is restricted (e.g. multiple non-cooperating 
providers in the same cable bundle), spectrum coordination 
(dynamic spectrum management) is applied ||3- Diagonal Pre¬ 
coding (DP) is a linear FEXT cancellation precoding adopted 
in VDSL2. It performs at the information theoretical limits in 
VDSL2 band where FEXT channel is much weaker than direct 
lines El. Non-linear FEXT cancellation based on Tomlinson- 
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Harashima Precoding (THP) ID provides significant gains over 
linear precoding in G.fast band where FEXT is often as strong 
as direct lines 0. 



This paper modifies the THP proposed in ||4l by introducing 
an ordering of successive interference pre-cancellation per¬ 
formed by CO that is optimized for non-cooperating CPEs. 
Ordering optimization has been already presented in 0 . 171 
for a different THP structure which we denote as Equal- 
Rate THP (ER-THP). ER-THP provides constant Signal-to- 
Noise Ratio (SNR) at each line. Any type of ordering can be 
concatenated with the proposed THP scheme. The ordering 
which maximizes the minimum rate on a single tone of 
Discrete Multi-Tone (DMT) modulation is V-BLAST (VB) 
ordering 0, ||9l. However, VB does not provide the maximal 
minimum of aggregated rates when applied on all DMT tones. 
We propose a novel Dynamic Ordering (DO) strategy that 
takes into account the asymmetry of G.fast channel statistics. 
The proposed scheme together with DO provides the maximal 
minimum rate of ~ 955 Mbps over tested 100 m long paper- 
insulated cable. The ordering arbitrary adjusting general user 
demands require considerable computation power IfTOl . com¬ 
plexity of DO is from this perspective negligible. 

Notation: Bold upper- and lower-case letters describe ma¬ 
trices and column vectors. [A],^ = denotes the ijth element 
of matrix A. Letters Z,Zj,R,C refer to integers, complex 
integers, real and complex numbers, respectively. We denote 
matrix inversion, transposition and conjugate transposition as 
(*)“', (a)^ , (*)^. Symbol = denotes equality from definition. 

Organization: Section [H] describes general THP scheme and 
its properties in which reference method 14) is defined. The 
proposed scheme is described in Sec. [Ill] and some ordering 
















Fig. 2. General Tomlinson-Harashima Precoding scheme. 


strategies are in Sec.|I3 Comparison to ER-THP, numerical 
results and conclusions are content of Sec. rVlIVllandlVlIl 

II. System Model 

A. Downstream Channel Model 

We assume centralized transmission from CO to non¬ 
cooperating CPEs in DS as shown in Eig.[T] DMT is employed 
to turn the frequency selective channel into a set of frequency 
flat orthogonal channels. On each tone, the signal received by 
L CPEs y e is modeled as 

y = Hx H-w, (1) 

where x = {x\,... ,xlY s denotes the transmitted signal 
vector and w e is the AWGN. The cable bundle is assumed 
to contain only the L lines. We avoid tone indexing to simplify 
the notation. Main diagonal elements /i„ = [H],, of channel 
matrix H e characterize insertion loss of direct lines, 

and off-diagonal elements hjj - [H]^ with i + j characterize 
EEXT. The channel is static and assumed to be known at the 
transmitter. 

B. General THP Scheme and Basic Properties 

We describe considered THP schemes in the common 
framework in Pig.|2] Linear block E represents the ordering (or 
later assumed lattice reduction). The feedback loop consisting 
of non-linear modulo Ft- block and linear block given by 
lower triangular matrix B with units along the main diagonal 
implements the inversion of B while reducing transmitted 
power by modulo F,-. F is a feedforward hlter which also 
ensures transmitted signal to satisfy energy constraints. Diag¬ 
onal matrix G describes linear operations performed by non¬ 
cooperating receivers. Let the input to the precoding chain be 
vector a = [ai,... ,aL]^ 6 forming data symbols and the 
output be decision variable vector y 6 C^. 

1) Linearized Scheme: Block Ft- is the modulo function 
over base t with origin shifted by Vz applied individually along 
each dimension of the input x. Particularly, 

FT[r] L(x + V2)modT - Vz, r 6 R, 

FtW = FT[5l{r)] H-jFjBlr)], x 6 C, 

^.[x] = [FT,[xl],...,FTJx2]]^ xeC^ 

Every modulo reminder equals to the input minus an integer 
multiple of base t such that the reminder is lower than t. 
Therefore Ft-[x] = x - d, where d is a vector such that x - 
d 6 [-^1/2, ■^i/z) X ■■■ X [-’■i/zAVz) ■ The ith component of d 
is di 6 T,Zj where Zj = Z -H jZ denotes complex integers and 
X = [ti, ..., is a vector of thresholds. Eigure [3 shows the 
linearized scheme where Ft- is replaced by additive term -d. 
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Fig. 3. Linearized general THP scheme. 
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Fig. 4. Threshold r for several considered square-shape QAM constellations. 

2) Zero-Forcing Condition: Zero-Forcing (ZF) precoding 
inverts the channel by eliminating the crosstalk such that 

y = GHFB 'Ea - GHFB 'd -h Gw - d (2) 

equals to input data a plus noise. The whole chain of linear 
blocks in Fig. [3] needs to fulfill ZF condition 

GHFB 'E = I, (3) 

with I being L x L identity matrix. Condition a implies 
GHFB-' = E-' which leads to y = a - E^'d + Gw - d. We 
obtain ZF property y = a -t- Gw when 

E-'d-Hd = 0 (4) 

which is realized by a proper design of thresholds r and r. 

3) Modulo Threshold t: Size of t is chosen to wrap 
constellations within r x t frame such that the distance from 
the edge point to the boundary is half of minimal distance 
dmia- Figure |4] shows the frame for several QAM constellations 
considered in this paper. It is straightforward to verify that 
square-shaped QAM constellations including odd-bit cardinal¬ 
ity variants (black points in Fig.|4|i have t = yfMd^i^i. Table |I] 
lists values of t considered here. 

4) Per-Line Power Constraint: Transmitted power on each 
line is constrained not to overcome a specified limit. Without 
loss of generality, we use constellations normalized to unit 
mean symbol energy E [|a,pj = 1 (£[★] denotes the statistical 
expectation) for which the energy limit implies 

£[W']<1. (5) 

We need to keep in mind and downscale constellations to pre¬ 
compensate energy increase AE due to modulo Ft-. Square- 
shape even-bit QAM constellations have AE ^ t^l{M-\) lU. 
The same formula holds for square-shape odd-bit cardinalities 
if twice higher value of M is used. For example, AE is the 
same for 2QAM and 4QAM as shown in Table D Odd-bit 
constellations with a square shape have signihcantly lower AE 
than popular cross-shaped constellations 0. 


TABLE I 

Threshold t and energy increase AF” due to modulo Ft- eor considered 

UNIT-MEAN SQUARE-SHAPE M-QAM CONSTELLATIONS. 


M 

2,4 

8,16 

32,64 

128,256 

512,1024 

2048,4096 

T 

2.83 

2.53 

2.47 

2.45 

2.45 

2.45 

AE [dB] 

1.25 

0.28 

0.068 

0.017 

0.0042 

0.0011 































C. Reference THP Scheme and Basic Properties 

The proposed THP scheme enhances the scheme described 
in m by ordered QR decomposition. Reference scheme H is 
described by definition of the blocks in Fig. [3 as 

E = I, B = diag(Rr*R", F = Q, G = diag(Rr\ (6) 

where diag(R)“' = diag ,... with diagonal compo¬ 
nents ra = [R]„ . Unitary matrix Q and upper-triangular R 
follow from QR decomposition of transposed channel matrix 

H" = QR. (7) 

Let us confirm basic properties introduced in Sec. lII-BI The 
reference scheme fulfills ZF condition (O 

GHFB^'E = diag(R)^'R"Q"QR^"diag(R)I = I. (8) 

Transmitted signal meets per-line energy constraint © 


B. Proposed TRIP Scheme and Basic Properties 

The ordered THP scheme is given by the following matrices 

E = P^,B = diag(Rr'R", F = Q, G = Pdiag(Rr'P^. (14) 

Key observation is that if diagonal matrix has permuted rows 
and columns by the same permutation (as G in (fT4]i ') then it 
remains diagonal and can be performed by non-cooperating 
receivers. Now, we show that ZF condition © is satisfied 

GHFB 'E = Pdiag(R) 'P^PR"Q"QR "diag(R)P^ = I, 

since P^ describes inverse permutation and so PP^ = I. As 
in the reference scheme in Sec. lll-CI feedforward matrix F is 
unitary and therefore transmitted signal meets per-line energy 
constraint ©. The decision variable y = a H- Pdiag(R) *P^w 
implies output SNR at the /th line to be 


L L 

7=1 

where qij - [Q],^ • We use the fact that x is approximately 
uncorrelated © and the energy increase due to Ft- has been 
pre-compensated E = 1. This last equality follows from 
that unitary Q has unit-length rows. The decision variable 


y = a H- diag(R) w 

implies that output SNR at the /th line is 

(10) 

yi — ybase ' Lj, 

(11) 


where jhase is the baseline input SNR. The main diagonal com¬ 
ponents can attain different values providing different 

SNR at the each line. In this case, different bit-loading per- 
line is required as well as different modulo threshold r, 4 tj 
in X = {ti, ... ,tlY, see Table |T] for actual values. There is 
no integer precoding operation as E = I and so selection of 
T - T fulfills modulo condition ©. 

III. Ordered Tomlinson-Harashima Precoding 
A. Ordered QR Decomposition 

The ordered THP scheme proposed here incorporates the 
ordered QR decomposition of transposed channel 


Ti ~ ybase ' (15) 

where ybase denotes baseline SNR and p, is the fth element 
of permutation output (fOl l. Similarly to the reference THP, 
different values of main diagonal components require 

different bit-loading with thresholds r = [ti ,..., where 
generally t, + tj. Vector of thresholds r needs to be chosen 
to fulfill condition © which means r = Px = [t^,, ..., Tpf\^. 

Remark 1. Any type of ordering P can be concatenated with 
the proposed ordered THP and it is a degree of freedom to be 
exploited. The reference scheme is obtained for the ordering 
P = I which means that optimized ordering can only improve 
the performance of the reference scheme. 

IV. Optimized Ordering of THP in G.fast Downstream 

There is a rich number of possible orderings (see IfTTl 
and references therein) to be concatenated with the scheme 
proposed here. Generally, different orderings lead to different 
SNR at each line O. Optimal selection is a multi-objective 
optimization problem where utility target considering fairness 
has significant impact on the result. We mainly focus on 
max-min fairness by maximizing the minimum rate and thus 
provide the same quality of service to each CPFs, although 
we discus sum-rate and simple combination of both as well. 


H" = QRP^ (12) 

into the reference scheme (Sec. lII-CT l. Permutation matrix P 
describes arbitrary permutation [ 1,...,L]^ —> {pi,... ,ppY as 



1 


Pi 


Cpi 

p 


= 


, where P = 



L 


PL 


. ^Pl - 


(13) 


and e, denotes a row vector with 1 in the /th position and 0 
elsewhere. Note that PX denotes permutation of rows of X 
and XP^ permutation of columns since XP^ = ■ 


A. V-BLAST (VB) Ordering 

The ordering strategy introduced in |[8l is the optimal 
max-min fair ordering maximizing the minimum SNR. The 
algorithm requires multiple calculations of channel matrix 
pseudo-inverse and so its complexity is much higher than 
the complexity of closely-related semi-optimal algorithm 0. 
Instead of ||8], we pragmatically use m since it performs close 
to the optimum without the computational burdens. 

VB ordering 0 is based on Gram-Schmidt (GS) QR 
decomposition of the transposed channel = H = QR. 
In the /th iteration, the algorithm choses the column vector 
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Fig. 5. Empirical probability density function of being weakest (i.e., being 
selected as the first to enter GS procedure in V-BLAST ordering). 
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where <h,q) = q«h denotes an inner product. The order in 
which h are chosen forms permutation matrix P in (fT2l l. This 
strategy (“weakest first”) leads to the ordering which maximize 
the minimum of elements and so SNR (fTSl l. 


Fig. 6. Let us assume a DMT system with 2 lines and 3 tones If/ll-i. 
Thick are strong lines with throughput of 6 bits and thin are weak lines with 
throughput of 1 bit. Entering GS procedure first is an advantage increasing 
the throughput by +1 bit and second decreasing the throughput by -1 bit. The 
minimum aggregated rate 2miii of V-BLAST (“weakest first”) is lower than 
DO (“aggregated minimum first”). Symbol 2, denotes aggregated rate at the 
(th line. We assume initial order at /i of DO to be given by VB. 


B. Inverse V-BLAST (IVB) Ordering 

IVB describes ordering with opposite approach than VB. In 
each GS iteration, always such a column vector h, is chosen 
for which diagonal element (fTdl l is maximal. It is a greedy 
maximization approach which maximize the sum-rate. We 
freely interchange sum-rate and mean-rate since the difference 
is just a scaling factor. The IVB ordering (“strongest first”) is 
also known as QR decomposition with pivoting lfT2ll . 
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Fig. 7. Frequency sharing between two extrema types of orderings DO 
(maximizing minimum) and IVB (maximizing sum-rate). Pai'ameter Bqo 
describes bandwidth assigned to DO. 


C. Dynamic Ordering (DO) 

Although VB is the optimal max-min fair ordering, we do 
not obtain equal rates when aggregated over multiple DMT 
tones in numerical results in Fig.[TT] So, if on a single tone, 
we selected instead of VB ordering a different ordering in favor 
of the line with the minimal aggregated rate, we would obtain 
the higher minimum. It means that VB ordering is optimal 
on a single tone, but it does not reach the global optimum if 
applied independently on each tones. The reason behind is that 
G.fast channel does not have the same statistical properties on 
each line. Some lines are more often the weakest lines (being 
selected first by VB) due to asymmetric physical arrangement 
of twisted-copper pairs within the cable bundle as confirmed 
by numerical evaluation in Fig. |5] 

We propose DO strategy taking into account this statis¬ 
tical asymmetry providing the highest minimum aggregated 
rate. DO ordering is inspired by VB approach which states 
that being taken first into GS is an advantage. Instead of 
VB “weakest first” approach applied independently on each 
tone, we propose to take first the line with so far minimum 
aggregated rate (“aggregated minimum first”). DO is ordering 
with memory deciding the order inductively in sequence. If 
the ordering on tone index 1 to / - 1 has been already chosen, 
then DO orders the lines on the /th tone as the order of 
bit-loading aggregated over tones from 1 to / - 1. Figure |6] 
shows an illustrative example explaining why DO provides 
higher aggregated minimum than VB. Complexity of DO is 


negligible, since the order is given by cumulative summation 
performed once at the beginning of transmission. The order is 
computed outside of QR algorithm and can be connected to 
whatever type of QR implementation, not only the one based 
on GS as in 0. 

D. Frequency-Sharing Between DO and IVB 

IVB ordering maximizes sum-rate on a single tone as well 
as when applied independently on multiple tones, unlike in 
max-min case of VB and DO ordering. We propose a simple 
frequency-sharing between two extrema types of ordering IVB 
and DO to adjust the fairness among CPFs. By frequency¬ 
sharing, we mean similar concept as time-sharing, but in 
the frequency domain. We propose to divide bandwidth on 
lower and upper parts where we expect different behavior 
(e.g., diagonal dominant property is present only on lower 
frequencies as shown in Fig.fTOli. DO ordering is allocated to 
lower frequencies in case DO-IVB and to higher frequencies 
in case IVB-DO as shown in Fig.|2] Numerical evaluation in 
Fig. [8] shows that DO-IVB sharing achieves better results. 

V. Comparison with Ordered Equal Rate THP (ER-THP) 

Ordering optimization has been introduced in 0, Q for 
THP structure which we denote as ER-THP (term centralized 
THP is also used ini). Label ’equal rate’ corresponds to the 
feature that ER-THP provides constant SNR. The proposed 
THP scheme and ER-THP have essentially different structure. 
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Fig. 8. Frequency sharing between DO and IVB enables to adjust trade ofi 
between mean and min-rate. For instance, both min-rate and mean-rate equa’ 
to ~ 950 Mbps for DO-IVB sharing with Sdo = 212 MHz, but if we permil 
a slight decrease of minimum, we could have mean-rate ~ 975 Mbps while 
having minimum still ~ 950 Mbps (here Bdo - 170 MHz). Similai'ly, wher 
mean-rate is priority, Bdo - 125 MHz increases min-rate from ~ 760 Mbps 
(Sdo - OMHz) to ~ 875 Mbps keeping the same mean-rate. 


A. Ordered ER-THP Scheme and Basic Properties 
Ordered ER-THP is defined by the following matrices 

E = P^,B = R"diag(Rr',F = -Qdiag(Rr', G = gl, (17) 

g 

where ordered QR decomposition (fT^ is used. Automatic gain 
control scaling g establishes power constrain 0 SO 

^ 1=1 

with labeling F = Qdiag(R)“' and |^f].. = fij. The constraint 
is fulfilled by the following scaling using (2, cx))-mixed norm 
ll*ll2,oc as 



Fig. 9. Direct line and crosstalk power characteristics of tested 100 m long 
paper-insulated cable. Values are smoothen by average over 1 MHz bin. 


Frequency / = 5 [MHz] 



Frequency / = 100 [MHz] 



Fig. 10. Diagonal dominant property of channel matrix |Hp [dB] at carrier 
frequency 5 [MHZ] is not present at higher frequency 100 [MHz]. 


{r„)^j) is again preferable. Reference IIT4l shows that even 
smaller value of g^ is obtained with LR QR decomposition 

H" = QRT^', (22) 


1=1 1=1 

Notice, that average power constrain assumed in 0, Q leads 
to the scaling with Frobenius norm ||★||;^ as g^ — i/l||f||^ = 
i/itr^FF^) . We confirm that ER-THP meets ZF condition 0 

GHFB 'E = gIPR"Q"iQdiag(R)-'diag(R)R = I. 

8 

The decision variable y = a -i- gw implies output SNR to be 


7i = ybase ■ V/. (20) 

where jbase denotes baseline SNR. Constant SNR yields 
the same bit-loading and the same modulo threshold r = 
[t, ..., t]^ on every line, therefore vector of thresholds t — t 
fulfills modulo condition (|4]i. 


B. VB Ordered and Lattice Reduced (LR) ER-THP 

Performance of ER-THP is given by scaling factor g^. Using 
inequality < i/min, r? , we rephrase (O as 


= max^ 
1=1 


\qi 


— < max ^ 
j ' J=i 


kvf 1 


mini “'"-t qk 


( 21 ) 


We see that g^ is minimized when min^ r^^ is as large as 
possible, therefore VB ordering (maximizing the minimum of 


where QR decomposes reduced channel as H = QR where 
reduced channel is H = H^T and T is a unimodular integer 
matrix. LR ER-THP is given by (fTTT i using decomposition (l22l i 
where E = T^. We use familiar LEE implementation of LR 
with moderate algorithm complexity parameter <5 = 3/4 US). 

Unfortunately, LR decomposition (l22l l cannot be used in 
the proposed THP scheme (fT4l l with E = T^, as matrix 
G = T“^ diag(R)“^T^ (unlike in the case of ordering) is 
not diagonal anymore and thus cannot be performed by non¬ 
cooperating CPEs. For the sake of comparison, we consider 
ER-THP scheme enhanced by both LR and VB ordering with 
complexity parameter set to an extreme value 5=1. The 
scheme has impractical implementation complexity but gives 
the highest SNR as confirmed by simulations in Fig.fTTI 

VI. Performance Evaluation over G.fast Channel 

A. Tested 100m Long Paper-Insulated Cable 

Figures l9] and [TOl show strong FEXT of 100m long paper- 
insulated G.fast cable where diagonal dominant property of 
VDSL2 is not present any more. 

B. Evaluation Procedure and Simulation Parameters 

Bit-loading at the /th line is computed according to Elby 
insertion of SNR y, (fTTTi . (fTSl) . (l20l) into the gap formula 

= Llog 2 (l+^'/r)J, for bie[2,l2] 


(23) 
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Fig. 11. Aggregated rates of considered precoding schemes over tested 
100m long paper-insulated cable. The acronyms and basic rate statistics of 
considered precoding schemes are summaiized in Table flUl 


TABLE II 

Simulation parameters (H 


Transmit PSD 

-76 [dBm/Hz] 

Coding Gain 

5[dB] 

Noise PSD 

-140 [dBm/Hz] 

Shannon Gap 

9.8 [dB] 

Band 

2.1 -212 [MHz] 

Bit Loading 

2 - 12 [bits] 

Tone Spacing A/ 

51.750 [kHz] 

Framing Overhead 

12 % 

Margin 

6 [dB] 




and bi - 0 otherwise, where symbol [★] denotes floor opera¬ 
tion and gap T - Shannon gap-tmargin-coding gain [dB]. Ag¬ 
gregated rates are obtained by summation of bit-loading (l2Bt 
over all DMT tones multiplied by A/ (1 - Framing Overhead), 
where A/ denotes the tone spacing. Considered parameters are 
listed in Table |II] We use bit-loading modification described 
in Is) to incorporate energy increases Afs (shown in Table |I]l 
due to modulo r,-. The algorithm allocates bits according to 
(l23b and then recomputes SNR corrected by AE and update 
bit-allocation accordingly. 

C. Numerical Results 

Numerical results in Fig. [TT] and Table Hill compare several 
FEXT cancellation methods in DS. Non-linear precoding 
based on THP or ER-THP clearly outperforms linear DP 
precoding used in VDSL2 [|3]. We confirm that ER-THP 
provides constant aggregated rates to all users, where gains 
by VB ordering and LR are significant. THP with un-equal 
rates provides higher sum-rate than ER-THP. As expected, VB 
ordered THP considerably increase min-rate and IVB ordered 
THP considerably increase sum-rate. Proposed DO ordered 
THP provides the highest aggregated minimum rate among 
all considered methods. The achieved rates are fairly stable 
vs. line index being close to G.fast target of 1 Gbps. 

VII. Conclusion 

Contribution of this paper is two fold, a new ordered THP 
scheme and a novel Dynamic Ordering (DO) strategy has 
been proposed, which together leads to the highest aggregated 
minimum rate in G.fast downstream. Unlike existing ordered 


TABLE III 

Mean and minimum aggregated rates of several FEXT cancellation 
PRECODING SCHEMES DEPICTED IN FiG. lllI ThE RATES ARE IN Mbits/s. 


Acronym 

Precoding Scheme 

[mean,min]-rate 

DP 

Diagonal Precoding 

[552,432] 

THP 

Tomlinson-Harashima Precoding 

[970,678] 

THP-VB 

THP using VB ordering 

[947,907] 

THP-IVB 

THP using Inverse VB ordering 

[990,763] 

THP-DO 

THP using Dynamic Ordering 

[956,955] 

ER-THP 

Equal-Rate THP 

[732,732] 

ER-THP-VB 

ER-THP using VB ordering 

[840, 840] 

ER-THP-LR 

ER-THP using Lattice Reduction 

[874, 874] 

ER-THP-LRVB 

ER-THP-LR using VB ordering 

[889, 889] 


Equal Rate (ER) THP scheme, the proposed scheme better 
adapts to asymmetric channel statistics of G.fast channel. 
Although the results are related to concrete G.fast settings, 
the proposed scheme has universal application in general 
multiple-input multiple-output systems including wireless sce¬ 
nario (e.g., paper lfT3]l shows that sum-rate of THP is always 
higher or equal than sum-rate of ER-THP when the same type 
of channel matrix decomposition is considered). 
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